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t~| , Abstract 

We use the method of thermal QCD sum rules to investigate the effects of temperature 
on the neutron electric dipole moment d n induced by the vacuum 0-angle. Then, we 
analyze and discuss the thermal behaviour of the ratio | % | in connection with the 
restoration of the CP-invariance at finite temperature. 
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1 Introduction 



The CP symmetry is, without doubt, one of the fundamental symmetries in nature. 
Its breaking still carries a cloud of mystery in particle physics and cosmology. Indeed, CP 
symmetry is intimately related to theories of interactions between elementary particles and 
represents a cornerstone in constructing grand unified and supersymmetric models. It is also 
necessary to explain the matter-antimatter asymmetry observed in universe. 

The first experimental evidence of CP violation was discovered in the K — K mixing and 
kaon decays |1[ . According to the CPT theorem, CP violation implies T violation. The latter 
is tested through the measurement of the neutron electric dipole moment (NEDM)d n . The 
upper experimental limit gives confidence that the NEDM can be another manifestation of 
CP breaking. To investigate the CP violation phenomenon many theoretical models were 
proposed. In the standard model of electroweak interactions, CP violation is parametrized by 
a single phase in the Cabbibo Kobayashi Maskawa (CKM) quark mixing matrix 0. Other 
models exhibiting a CP violation are given by extensions of the standard model; among 
them, the minimal supersymmetric standard model MSSM includes in general soft complex 
parameters which provide new additional sources of CP violation |3|, |4[]. 

CP violation can be also investigated in the strong interactions context through QCD 
framework. In fact, the QCD effective lagrangian contains an additional CP-odd four dimen- 
sional operator embedded in the following topological term: 

U = O^G^GT, (1.1) 

where G^ u is the gluonic field strength, G^ u is its dual and a s is the strong coupling constant. 
The G^G^ quantity is a total derivative, consequently it can contribute to the physical 
observables only through non perturbative effects. The NEDM is related to the #-angle by 
the following relation : 

e ,m 2.7 xlO- 16 ^ [5] 

dn WMr) ~ { 5.2x 10-^ i (L2) 

and consequently, according to the experimental measurements d n < 1.1 x 10~ 25 ecm 0, the 
9 parameter must be less than 2 x 10~ 10 ||. The well known strong CP problem consists 
in explaining the smallness of 6. In this regard, several scenarios were suggested. The most 
known one was proposed by Peccei and Quinn [1(J and consists in implementing an extra Ua(1) 
symmetry which permits a dynamical suppression of the undesired #-term. This is possible 
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due to the fact that the axial current J§ is related to the gluonic field strength through the 
following relation Ji* = ^G^ V G^ V ' . The breakdown of the J7a(1) symmetry gives arise to 
a very light pseudogoldstone boson called axion. This particle may well be important to the 
puzzle of dark matter and might constitute the missing mass of the universe ]TT] . 

Motivated by: (a) the direct relation between the Wangle and NEDM d n , as it was demon- 
strated firstly in || via the chiral perturbation theory and recently in [|T2"| , |13|| within QCD 
sum rules formalism; (b) the possibility to restore some broken symmetries by increasing the 
temperature; we shall use the QCD sum rules at T 7^ to derive thermal dependence of 
the ratio | % |. Then we study its thermal behaviour at low temperatures and discuss the 
consequences of temperature effects on the restoration of the broken CP symmetry. 

This paper is organized as follows: Section 2 is devoted to the calculations of the NEDM 
induced by the 9 parameter from QCD sum rules. In section 3, we show how one intro- 
duces temperature in QCD sum rules calculations. We end this paper with a discussion and 
qualitative analysis of the thermal effects on the CP symmetry. 



2 NEDM from QCD sum rules 



In the two later decades, QCD sum rules a la SVZ [|1J] were applied successfully to the 
investigation of hadronic properties at low energies. In order to derive the NEDM through 
this approach, many calculations were performed in the literature |TR ESI . One of them, which 



turns out to be more practical for our study, has been obtained recently in [12, 13]. It consists 



in considering a lagrangian containing the following P and CP violating operators: 

Lp,cp = -6 q m, £ q f i l5 q f + 9^-G, u G^. (2.1) 

8 q and 9 are respectively two angles coming from the chiral and the topological terms and 
m* is the quark reduced mass given by m * = m U +m d • ^ ne au thors of [TBj start from the two 
points correlation function in QCD background with a nonvanishing 9 and in the presence of 
a constant external electomagnetic field F^ v : 



II(g 2 ) = % J d 4 xe iqx < 0\T{rj{x)fj{0)}\0 > e , F • (2.2) 



rj(x) is the interpolating current which in the case of the neutron reads as ||T6f : 

V = 2e ahc {{d T a C lh u b )d c + (3(d^Cu b ) 75 d c }, (2.3) 
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where (3 is a mixing parameter. Using the operator product expansion (OPE), they have first 
performed the calculation of II (g 2 ) as a function of matrix elements and Wilson coefficients and 
then have confronted the QCD expression of II (g 2 ) to its phenomenological parametrisation. 
Il(g 2 ) can be expanded in terms of the electromagnetic charge as| 



II(g 2 ) =U^\q 2 ) + eU^(q\F^) + 0(e 2 ). (2.4) 

The first term nJ°)(g 2 ) is the nucleon propagator which include only the CP-even parameters 
19 , |1~8|, while the second term Il^^g 2 , F^ u ) is the polarization tensor which may be expanded 



through Wilson OPE as: C n < 0|gPg|0 >e,F, where T is an arbitrary Lorentz structure 
and C n are the Wilson coefficient functions calculable in perturbation theory ||19|| . From this 



expansion, one keeps only the CP-odd contribution piece. By considering the anomalous axial 



current, one obtains the following 9 dependence of < 0|gPg|0 >q matrix elements |T2 



m q < 0|gTg|0 > = im*6 < 0|gTg|0 >, (2.5) 

where m q and m* are respectively the quark and reduced masses. The electromagnetic depen- 
dence of these matrix elements can be parametrized through the implementation of the k, x 



and £ susceptibilities defined as |18| : 



< 0\qa^q\0 > F = yF^ < 0|gg|0 > 

g <0\qG^q\0 > F = kF^ u <0\qq\0 > (2.6) 
2g < 0|gG^g|0 > F = ^F" v < 0|gg|0 > . 



Putting altogether the above ingredients and after a straightforward calculation ||13| , the 
following expression of Il^^g 2 , F^) for the neutron is derived: 

n ("9 2 ) = -|^<0|gg|0>{F ( r,g}[x(/?+l) 2 (4e,-e u )ln( A " 



-409 - l) 2 e d (l + i(2« + 0)(ln(^) - 1) ' 



2 - 



4 V ~" ' *" y ~ s (i 2 IR ' -g 2 

|((4/3 2 - 4/3 + 2)e d + (3{3 2 + 2(3+ 1)^) J-...], (2.7) 
2 — q A 



where 9 = 9 + 9 q is the physical phase and g = g M 7 At . 

The QCD expression (2.7) will be confronted to the phenomenological parametrisation U Phen (—q 2 ) 
written in terms of the Neutron hadronic properties. The latter is given by: 

n-(- 9 >) = { ^ Kl ^ + ^ + ...), (2.8, 
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where m n is the neutron mass, e q is the quark charge. A and A 2 , which originate from the 
phenomenological side of the sum rule, represent respectively a constant of dimension 2 and 
the neutron coupling constant to the interpolating current r)(x). This coupling is defined via 
a spinor v as < 0\r](x)\n >= Xv. 



3 QCD sum rules at finite temperature 

The introduction of finite temperature effects may provide more precision to the phe- 
nomenological values of hadronic observables. Within the framework of QCD sum rules, the 
T-evolution of the correlation functions appear as a thermal average of the local operators in 
the Wilson expansion pU|, ET0, UM. Hence, at nonzero temperature and in the approximation of 



the non interacting gas of bosons (pions), the vacuum condensates can be written as : 

< & > T =< 0*>+J < n(p)\OMp) > (3-1) 



where e = yp 2 + m 2 , Ub = ^r[^ s ^ ne Bose-Einstein distribution and < % > is the standard 
vacuum condensate (i.e. at T=0). In the low temperature region, the effects of heavier 
resonances (r = K, rj, ..etc) can be neglected due to their distibution functions ~ e t 
To compute the pion matrix elements, we apply the soft pion theorem given by: 

1 • 771 2 

< n(pW\n(p) >= -- < 0\[F^ [F 5 a ,OW > +0(^f), (3.2) 

where A is a hadron scale and F£ is the isovector axial charge: 

F 5 a = J d 3 xq(x) l0l5 ^q(x). (3.3) 

Direct application of the above formula to the quark and gluon condensates shows that [ 22 , : 
(i) Only < qq > is sensitive to temperature. Its behaviour at finite T is given by: 

< qq > T ~ (1 - ^P-) < qq >, (3.4) 



where y?(T) = j^Bi^f-) with B(z) = Jy dy^^— and /„- is the pion decay constant 
(/,,. ~ 93MeV). The variation with temperature of the quark condensate < qq >t results in 
two different asymptotic behaviours, namely: 

< qq > T ~ (1 - ^) < qq> for < 1, and < qq > T ^ (1 - f^e^) < qq > for ^ > 1. 
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(ii) The gluon condensate is nearly constant at low temperature and a T dependence occurs 
only at order T 8 . 

As usual, the determination of the ratio ~ sum rules at non zero temperature is now easily 
performed through two steps. In the first step, we apply Borel operator to both expressions 
of the Neutron correlation function shown in Eqs. (2.7) and (2.8), where finite temperature 
effects were introduced as discussed above. Next step, by invoking the quark-hadron duality 
principle, we deduce the following relation of the 6 induced NEDM: 

T< T > = -^AflrlW)* 1 -^' < ^ > We.-e^J^e^e^, (3.5) 

where M represents the Borel parameter. Note that we have neglected the single pole contri- 
bution entering via the constant A, as suggested in 

The expression (3.5) is derived with j3 — 1 which is more appropriate for us since it suppresses 
the infrared divergences. In fact, the Ioffe choice (/3 = —1) which is rather more useful for the 
CP even case, removes the leading order contribution in the sum rules (2.7). The coupling 
constant A^(T) and the neutron mass M n (T) which appear in (3.5) were determined from the 
thermal QCD sum rules. For the former, we consider the q sum rules in |2"5| , [^J with (3 = 1 
and then we extract the following explicit expression of A^(T): 

Xl(T) = {— ^M 6 + — ^M 2 < —G 2 >}{1 - (1 + %^)^} e 4^ (3.6) 

nK 1 X 8(2^) 4 16(2tt) 2 tv 11 y Ml ! 16 1 y J 



Within the pion gas approximation, Eletsky has demonstrated in [gfj that inclusion of the 
contribution coming from the pion-nucleon scattering in the nucleon sum rules is mandatory. 
The latter enters Eq.(3.6) through the coupling constant g W NN, whose values lie within the 



range 13.5-14.3 |27 



Numerical analysis is performed with the following input parameters: the Borel mass has 
been chosen within the values M 2 = 0.55 — 0.7GeV 2 which correspond to the optimal range 
(Borel window) in the % sum rule at T = |13|]. For the \ an d £ susceptibilities we take 



X = — 5.7 ± 0. 6GeV 2 ]28j and £ = —0.74 ±0.2 As to the vacuum condensates appearing 



in (3.5), we fix < qq > and < G > to their standard values |14 . 

4 Discussion and Conclusion 



In the two above sections, we have established the relation between the NEDM and 6 
angle at non zero temperature from QCD sum rules. Since the ratio %■ is expressed in terms 
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of the pion parameters /„., and of g^NN, we briefly recall the main features of their thermal 
behaviour. Various studies performed either within the framework of the chiral perturbation 
theory and/or QCD sum rules at low temperature have shown the following features: 
(i) The existence of QCD phase transition temperature T c which signals both QCD deconfine- 



ment and chiral symmetry restoration |2(J, |3T 



(ii) f n and g n NN have very small variation with temperature up to T c . So, we shall assume 
them as constants below T c . However, they vanish if the temperature passes through the 



critical value T c pO 



(iii) The thermal mass shift of the neutron and the pion is absent at order 0(T 2 ) |23], [2"T| ]. 
5M n shows up only at the next order T 4 , but its value is negligible [26 ]. 
By taking into account the above properties, we plot the ratio defined in Eq. (3.5) as a func- 
tion of T. From the figure, we learn that the ratio | % | survives at finite temperature and 
it decreases smoothly with T (about 16% variation for temperature values up to 200 MeV). 
This means that either the NEDM value decreases or 9 increases. Consequently, for a fixed 
value of 9 the NEDM decreases but it does not exhibit any critical behaviour. Furthermore, 
if we start from a non vanishing 9 value at T = 0, it is not possible to remove it at finite 
temperature. We also note that | | grows as M 2 or \ susceptibility increases. It also 
grows with quark condensate rising. However this ratio is insensitive to both the £ suscep- 
tibility and the coupling constant g W NN- We notice that for higher temperatures, the curve 
| ^& |= /(|r) exhibits a brutal increase justified by the fact that for temperatures beyond the 
critical value T c , at which the chiral symmetry is restored, the constants f n and g^NN become 
zero and consequently from Eq(3.5) the ratio % behaves as a non vanishing constant. The 
large difference between the values of the ratio for T < T c and T > T c maybe a consequence 
of the fact that we have neglected other contributions to the the spectral function, like the 
scattering process TV + n — > A. These contributions , which are of the order T 4 , are negligible 
in the low temperture region but become substantial for T > T c . Moreover, this difference 
may also originate from the use of soft pion approximation which is valid essentially for low T 
(T < T c ). Therefore it is clear from this qualitative analysis, which is based on the soft pion 
approximation, that the temperature does not play a fundamental role in the suppression of 
the undesired #-term and hence the broken CP symmetry is not restored, as expected. This is 
not strange, in fact it was shown that more heat does not imply automatically more symmetry 



33], |36|. Moreover, some exact symmetries can be broken by increasing temperature [|32], pqj . 
The symmetry non restoration phenomenon, which means that a broken symmetry at T=0 
remains broken even at high temperature, is essential for discrete symmetries, CP symmetry 
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in particular. Indeed, the symmetry non restoration allows us to avoid wall domains inherited 
after the phase transition and to explain the baryogenesis phenomenon in cosmology [j35 ]. 



Furthermore, it can be very useful for solving the monopole problem in grand unified theories 
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Figure Captions 

Figure: Temperature dependence of the ratio 
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